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§1. Introduction 

This paper completes the proof of the celebrated Boltzmann-Sinai Ergodic Hy- 
pothesis. In a loose form, as attributed to L. Boltzmann back in the 1880's, this 
hypothesis asserts that gases of hard balls are ergodic. In a precise form, which is 
due to Ya. G. Sinai in 1963 [Sin(1963)], it states that the gas of TV > 2 identical 
hard balls (of "not too big" radius) on a torus T'^ , u > 2, (a //-dimensional box with 
periodic boundary conditions) is ergodic, provided that certain necessary reductions 
have been made. The latter means that one fixes the total energy, sets the total 
momentum to zero, and restricts the center of mass to a certain discrete lattice 
within the torus. The assumption of a not too big radius is necessary to have the 
interior of the configuration space connected. 

Sinai himself pioneered rigorous mathematical studies of hard ball gases by prov- 
ing the hyperbolicity and ergodicity for the case N = 2 and u = 2 in his seminal 
paper [Sin(1970)], where he laid down the foundations of the modern theory of 
chaotic billiards. Then Chernov and Sinai extended this result to {N = 2, z/ > 2), 
as well as proved a general theorem on "local" ergodicity applicable to systems of 
N > 2 balls [S-Ch(1987)]; the latter became instrumental in the subsequent stud- 
ies. The case N > 2 is substantially more difficult than that of N = 2 because, 
while the system of two balls reduces to a billiard with strictly convex (spherical) 
boundary, which guarantees strong hyperbolicity, the gases of > 2 balls reduce to 
billiards with convex, but not strictly convex, boundary (the latter is a finite union 
of cylinders) - those are characterized by very weak hyperbolicity. 

Further development has been mostly due to A. Kramli, D. Szasz, and the present 
author. We proved hyperbolicity and ergodicity for N = 3 balls in any dimension 
[K-S-Sz(1991)] by exploiting the "local" ergodic theorem of Chernov and Sinai [S- 
Ch(1987)], and carefully analyzing all possible degeneracies in the dynamics to 
obtain "global" ergodicity. We extended our results to = 4 balls in dimension u > 
3 next year, and then I proved the ergodicity whenever N <v (this covers systems 
with an arbitrary number of balls, but only in spaces of high enough dimension, 
which is a restrictive condition). At this point, the existing methods could no 
longer handle any new cases, because the analysis of the degeneracies became overly 
complicated. It was clear that further progress should involve novel ideas. 

A breakthrough was made by Szasz and myself, when we used the methods of 
algebraic geometry [S-Sz(1999)]. We assumed that the balls had arbitrary masses 
mi, . . . ,mN (but the same radius r). Now by taking the limit m^v — 0, we were 
able to reduce the dynamics of balls to the motion of iV — 1 balls, thus utilizing a 
natural induction on N. Then algebro- geometric methods allowed us to effectively 
analyze all possible degeneracies, but only for typical (generic) vectors of "external" 
parameters (mi, . . . ,mN, r); the latter needed to avoid some exceptional submani- 
folds of codimension one, which remained unknown. This approach led to a proof of 
full hyperbolicity (but not yet ergodicity) for all A'" > 2 and ly >2, and for generic 
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(mi, . . . ,mN,r), see [S-Sz(1999)]. Later I simplified the arguments and made them 
more "dynamical" , which allowed me to obtain full hyperbolicity for hard balls with 
any set of external geometric parameters (mi, . . . , mAr, r) [Sim(2002)]. Thus, the 
hyperbolicity has been fully established for all systems of hard balls on tori. 

To upgrade the full hyperbolicity to ergodicity one needs to refine the analy- 
sis of the aforementioned degeneracies. For hyperbolicity, it was enough that the 
degeneracies made a subset of codimension > 1 in the phase space. For ergodic- 
ity, one has to show that its codimension is > 2, or to find some other ways to 
prove that the (possibly) arising codimension-one manifolds of non- sufficiency are 
incapable of separating distinct ergodic components. The latter approach will be 
pursued in this paper. In the paper [Sim(2003)] I took the first step in the direction 
of proving that the codimension of exceptional manifolds is at least two: I proved 
that the systems of > 2 balls on a 2D torus (i.e., u = 2) are ergodic for typical 
(generic) vectors of external parameters (mi, . . . ^rriN^r). The proof again involves 
some algebro- geometric techniques, thus the result is restricted to generic parame- 
ters (mi, . . . , mAr; r). But there was a good reason to believe that systems in u > S 
dimensions would be somewhat easier to handle, at least that was indeed the case 
in early studies. 

Finally, in my recent paper [Sim(2004)] I was able to further improve the algebro- 
geometric methods of [S-Sz(1999)], and proved that for any N > 2, u > 2 and for 
almost every selection (mi, . . . jUIn; r) of the external geometric parameters the 
corresponding system of N hard balls on T'^ is (fully hyperbolic and) ergodic. 

In this paper I will prove the following result. 

Theorem. For any integer values N > 2, f > 2, and for every (A^ + l)-tuple 
(mi, . . . , m^v, r) of the external geometric parameters the standard hard ball system 



Remark 1.1. The novelty of the theorem (as compared to the result in [Sim(2004)]) 
is that it applies to each (N + l)-tuple of external parameters (provided that the 
interior of the phase space is connected), without an exceptional zero-measure set. 

Remark 1.2. The present result speaks about exactly the same models as the 
result of [Sim(2002)], but the assertion of this new theorem is obviously stronger 
than that of the theorem in [Sim(2002)]: It has been known for a long time that, for 
the family of semi-dispersive billiards, ergodicity cannot be obtained without also 
proving full hyperbolicity. 

Remark 1.3. As it follows from the results of [C-H(1996)] and [0-W(1998)], all 
standard hard ball systems (M, {5'*}t£]R, //) (the models covered by the theorem) 
are not only ergodic, but they enjoy the Bernoulli mixing property, as well. 




is (fully hyperbolic and) ergodic. 
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The Organization of the Paper. In the subsequent section we overview the 
necessary technical prerequisites of the proof, along with the needed references to 
the literature. The fundamental objects of this paper are the so called "exceptional 
J-manifolds" : they are codimension-one submanifolds of the phase space that are 
separating distinct, open ergodic components of the billiard flow. In §3 we obtain 
the necessary linear estimations for the contraction coefficients of some specific 
tangent vectors in small, tubular neighborhoods of exceptional manifolds. By using 
the results of §3, in §4 we prove that at least one phase point of an exceptional 
J-manifold is actually sufficient (Main Lemma 4.5). In §5 we describe a fascinating, 
new, relaxed version of the Ansatz, sufficient for the proof of the Theorem on Local 
Ergodicity. In the same section we also prove this relaxed Ansatz by using the 
assumed full hyperbolicity of the semi-dispersive billiard system. 

Finally, in the closing section we complete the inductive proof of ergodicity (with 
respect to the number of balls TV) by utilizing Main Lemma 4.5 and the results of §5. 
Actually, a consequence of Main Lemma 4.5 will be that exceptional J-manifolds do 
not exist, and this will imply the the fact that no distinct, open ergodic components 
can coexist. 

Finally, a short appendix of this paper serves the purpose of making the reading 
of the proof easier, by providing a chart of the hierarchy of the selection of several 
constants playing a role in the proof of Main Lemma 4.5. 



§2. Prerequisites 



Consider the z/-dimensional {u > 2), standard, flat torus T*^ = W jll' as the 

vessel containing (> 2) hard balls (spheres) Bi, . . . , with positive masses 

mi, . . . jUiN and (just for simplicity) common radius r > 0. We always assume that 

the radius r > is not too big, so that even the interior of the arising configuration 

space Q (or, equivalently, the phase space) is connected. Denote the center of the 

ball Bi by qi e T'^, and let Vi = qi be the velocity of the i-th particle. We investigate 

the uniform motion of the balls Si, ... , S^v inside the container T'^ with half a unit 

1 1 

of total kinetic energy: E = - Xli=i "^^ll'^jll^ — 2' assume that the collisions 
between balls are perfectly elastic. Since — beside the kinetic energy E — the total 
momentum / = Xli^i '^i'^i ^ is also a trivial first integral of the motion, we 
make the standard reduction / = 0. Due to the apparent translation invariance of 
the arising dynamical system, we factorize the configuration space with respect to 
uniform spatial translations as follows: (q*!, . . . , q^) ~ {q\ -\- a, . . . ^q^ -\- a) for all 
translation vectors a G T^. The configuration space Q of the arising flow is then 

the factor torus ^(T*^)^ / = Y'^i^-'^) minus the cylinders 

Ci,j = {(gi, . . . , giv) e T'^(^-i): dist(g,, g,) < 2r} 
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{1 < i < j < N) corresponding to the forbidden overlap between the i-th and j-th 
spheres. Then it is easy to see that the compound configuration point 

9=(9i,...,giv)eQ = T-^^-'^ U C.j 

l<i<j<N 

moves in Q uniformly with unit speed and bounces back from the boundaries dCij 
of the cylinders Cij according to the classical law of geometric optics: the angle of 
reflection equals the angle of incidence. More precisely: the post-collision velocity 
v'^ can be obtained from the pre-coUision velocity v~ by the orthogonal reflection 
across the tangent hyperplane of the boundary dQ, at the point of collision. Here we 
must emphasize that the phrase "orthogonal" should be understood with respect to 
the natural Riemannian metric (the kinetic energy) = S^Li^ilM^ilP the 

configuration space Q. For the normalized Liouville measure fx of the arising flow 
{S^} wc obviously have d/j, = const ■ dq ■ dv, where dq is the Riemannian volume 
in Q induced by the above metric, and dv is the surface measure (determined by 
the restriction of the Riemannian metric above) on the unit sphere of compound 
velocities 



{N N ^ 

ivi,...,VN) eC^^f : E™i^* = Oand ^ ™»l 1^*1 1^ = 1 f • 
i=l i=l ) 

The phase space M of the flow {5*} is the unit tangent bundle Q x S'^"-'^ of the 
configuration space Q. (Wc will always use the shorthand notation d = v{N — 1) 
for the dimension of the billiard table Q.) We must, however, note here that at 
the boundary of Q one has to glue together the pre-coUision and post-collision 
velocities in order to form the phase space M, so M is equal to the unit tangent 
bundle Q x S*^"^ modulo this identification. 

A bit more detailed definition of hard ball systems with arbitrary masses, as well 
as their role in the family of cylindric billiards, can be found in §4 of [S-Sz(2000)] 
and in §1 of [S-Sz(1999)]. We denote the arising fiow by (M, {5*}teR, y). 

In the series of articles [K-S-Sz(1989)], [K-S-Sz(1991)], [K-S-Sz(1992)], [Sim(1992- 
I)], and [Sim(1992-II)] the authors developed a powerful, three-step strategy for 
proving the (hyperbolic) ergodicity of hard ball systems. First of all, these proofs 
are inductions on the number of balls involved in the problem. Secondly, the 
induction step itself consists of the following three major steps: 

Step I. To prove that every non-singular (i. e. smooth) trajectory segment S^°''^^Xq 
with a "combinatorially rich" (in a well defined sense) symbolic collision sequence is 
automatically sufficient (or, in other words, "geometrically hyperbolic", see below 
in this section), provided that the phase point xq does not belong to a countable 
union J of smooth sub-manifolds with codimension at least two. (Containing the 
exceptional phase points.) 
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The exceptional set J featuring this result is negligible in our dynamical consid- 
erations — it is a so called slim set. For the basic properties of slim sets, see again 
below in this section. 

Step II. Assume the induction hypothesis, i. e. that all hard ball systems with N' 
balls {2 < N' < N) are (hyperbolic and) ergodic. Prove that there exists a slim set 
5* C M with the following property: For every phase point xq G M \ the entire 
trajectory S^xq contains at most one singularity and its symbolic collision sequence 
is combinatorially rich, just as required by the result of Step I. 

Step III. By using again the induction hypothesis, prove that almost every singular 
trajectory is sufficient in the time interval {to, +oo), where to is the time moment 
of the singular reflection. (Here the phrase "almost every" refers to the volume 
defined by the induced Riemannian metric on the singularity manifolds.) 

We note here that the almost sure sufficiency of the singular trajectories (fea- 
turing Step III) is an essential condition for the proof of the celebrated Theorem 
on Local Ergodicity for semi-dispersive billiards proved by Chernov and Sinai [S- 
Ch(1987)]. Under this assumption, the result of Chernov and Sinai states that in 
any semi-dispersive billiard system a suitable, open neighborhood Uq of any suffi- 
cient phase point G M (with at most one singularity on its trajectory) belongs 
to a single ergodic component of the billiard flow (M, {5'*}teR, //). 

A few years ago Balint, Chernov, Szasz, and Toth [B-Ch-Sz-T(2002)] discovered 
that, in addition, the algebraic nature of the scatterers needs to be assumed, in 
order for the proof of this result to work. Fortunately, systems of hard balls are, by 
nature, automatically algebraic. 

In an inductive proof of ergodicity, steps I and II together ensure that there 
exists an arc-wise connected set C C M with full measure, such that every phase 
point a^o G C* is sufficient with at most one singularity on its trajectory. Then the 
cited Theorem on Local Ergodicity (now taking advantage of the result of Step III) 
states that for every phase point xq G C an open neighborhood Uq of xq belongs 
to one ergodic component of the flow. Finally, the connectedness of the set C and 
//(M \ C) = imply that the flow (M, {S'^jteM, /x) (now with N balls) is indeed 
ergodic, and actually fully hyperbolic, as well. 

The generator subspace Aij C {1 < i < j < N) oi the cylinder Cij 

(describing the collisions between the i-th. and j-th. balls) is given by the equation 

(2.1) Aij = [{q^, ...,qN)e {Wf : Qi = Qj} , 

see (4.3) in [S-Sz(2000)]. Its ortho-complement Lij C is then deflned by the 
equation 

(2.2) Lij = I (gi, . . . , gAr) G (R")^ : gfe = for 7^ i,j, and ruiqi + mjqj = [ , 
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see (4.4) in [S-Sz(2000)]. Easy calculation shows that the cylinder Ci^j (describing 
the overlap of the z-th and j-th balls) is indeed spherical and the radius of its base 

sphere is equal to r^j- = 2r.^^^^^, see §4, especially formula (4.6) in [S-Sz(2000)]. 

The structure lattice C C is clearly the lattice C = (Z^)^ = Z^''. 

Due to the presence of an additional invariant quantity / = X^i^i "^i'^^i; one 
usually makes the reduction J^^^ = and, correspondingly, factorizes the 
configuration space with respect to uniform spatial translations: 

(2.3) {qi,...,qN) ^ {qi + a,...,qN + a), aeT". 

The natural, common tangent space of this reduced configuration space is 



AT 



i=l 



(2.4) Z=^{vi,...,VN)e (R^f : J2 ^i^i = r = I n I = (-^)^ 
supplied with the inner product 




AT 



{v, v') = ^mi(wj, v[), 



see also (4.1) and (4.2) in [S-Sz(2000)]. 

Collision graphs. Let S't^'^^x be a nonsingular, finite trajectory segment with 
the collisions cri, . . . ,o"„ listed in time order. (Each is an unordered pair (z,j) 
of different labels i,j G {1, 2, . . . , A'"}.) The graph Q — (V, £) with vertex set 
V = {1,2,..., N} and set of edges £ = {ui, . . . , cr„} is called the collision graph 
of the orbit segment S^'^'^^x. For a given positive number C, the collision graph 
Q — (VjS) of the orbit segment S^'^'^^x will be called C-rich if Q contains at least 
C connected, consecutive (i. e. following one after the other in time, according to 
the time-ordering given by the trajectory segment S^'^'^^x) subgraphs. 

Trajectory Branches. We are going to briefly describe the discontinuity of the 

flow {S*} caused by a multiple collisions at time to- Assume flrst that the prc- 
coUision velocities of the particles are given. What can we say about the possible 
post-collision velocities? Let us perturb the pre-coUision phase point (at time 
to — 0) inflnitesimally, so that the collisions at ~ to occur at inflnitesimally different 
moments. By applying the collision laws to the arising flnite sequence of collisions, 
we see that the post-collision velocities are fully determined by the time-ordered list 
of the arising collisions. Therefore, the collection of all possible time-ordered lists of 
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these collisions gives rise to a finite family of continuations of the trajectory beyond 
to. They are called the trajectory branches. It is quite clear that similar statements 
can be said regarding the evolution of a trajectory through a multiple collision in 
reverse time. Furthermore, it is also obvious that for any given phase point G M 
there are two, a;-high trees and 71 such that (71) describes all the possible 
continuations of the positive (negative) trajectory S^^'°°'^xo {S^~°°'^^xo). (For the 
definitions of trees and for some of their applications to billiards, cf. the beginning 
of §5 in [K-S-Sz(1992)].) It is also clear that all possible continuations (branches) 
of the whole trajectory S^~°°'°°'^Xo can be uniquely described by all pairs (S_, B^) 
of infinite branches of the trees 71 and (i?_ C 71, G T^). 

Finally, we note that the trajectory of the phase point xq has exactly two 
branches, provided that S^xq hits a singularity for a single value t = to, and the 
phase point S^°xo does not lie on the intersection of more than one singularity 
manifolds. In this case we say that the trajectory of xo has a "simple singularity". 

Neutral Subspaces, Advance, and Sufflciency. Consider a nonsingular tra- 
jectory segment S't^'^lx. Suppose that a and b are not moments of collision. 

Definition 2.5. The neutral space N'o{S^"''^^x) of the trajectory segment S^"''^^x at 
time zero (a < < b) is defined by the following formula: 

AfoiS^'^'^^x) = {W eZ: 3{S > 0) such that Va e (-6,6) 

V {S"" {Q{x) + aW, V{x))) = ViS^'x) and V {S^ {Q{x) + aW, V{x))) = V{S^x)}. 

[Z is the common tangent space T^Q of the parallelizable manifold Q at any of its 
points while V{x) is the velocity component of the phase point x = {Q{x), V{x)).) 

It is known (see (3) in §3 of [S-Ch (1987)]) that AfoiS^'^'^^x) is a hnear subspace 
of Z indeed, and V{x) G N'o{S^'^'^^x). The neutral space N't{S^°''^^x) of the segment 
S^"''^^x at time t e [a, b] is defined as follows: 

Mt{S^'''b^x) =AAo (^S^''-''^-*^iS'x)') . 

It is clear that the neutral space can be canonically identified with 

A/o(5't"'^]x) by the usual identification of the tangent spaces of Q along the trajectory 
^(-oo,oo)^ (see, for instance, §2 of [K-S-Sz(1990)-I]). 

Our next definition is that of the advance. Consider a non-singular orbit segment 
S^°''^^x with the symbolic collision sequence S = (ci, . . . , c^), meaning that 5 X 
has exactly n collisions with dQ, and the z-th collision {1 < i < n) takes place at the 
boundary of the cylinder Co-.. For x — {Q,V) e M and W E Z, \\W\\ sufficiently 
smaU, denote Tw{Q,V) := {Q + W, V). 
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Definition 2.6. For any 1 < k < n and t e [a, b], the advance 

ak = a{ak): AftiS^'''^) ^ R 

of the collision ak is the unique linear extension of the linear functional au = 
a{ak) defined in a sufficiently small neighborhood of the origin of J\ft{S^°''^^x) in the 
following way: 

a{<Jk){W) := tk{x) - tk{S-'TwS'x). 

Here tk = tk{x) is the time of the A;-th coUision Ufe on the trajectory of x after 
time t = a. The above formula and the notion of the advance functional 

has two important features: 

(i) If the spatial translation {Q,V) ^ {Q + W, V) (W e M (S^'^'^^x)) is carried 
out at time t, then tk changes linearly in and it takes place just ak{W) units of 
time earlier. (This is why it is called "advance".) 

(ii) If the considered reference time t is somewhere between tk-i and tk, then 
the neutrality of W with respect to ak precisely means that 

W-ak{W)-V{x) e A,,, 

i. e. a neutral (with respect to the collision ak) spatial translation W with the 
advance ak{W) = means that the vector W belongs to the generator space A^^. 
of the cylinder (7^^ . 

It is now time to bring up the basic notion of sufficiency (or, sometimes it is 
also called geometric hyperbolicity) of a trajectory (segment). This is the utmost 
important necessary condition for the proof of the Theorem on Local Ergodicity for 
semi-dispersive billiards, [S-Ch(1987)]. 

Definition 2.7. 

(1) The nonsingular trajectory segment S^°''''^x (a and b are supposed not to be 
moments of collision) is said to be sufficient if and only if the dimension of 
A/;(5['*'*la;) (t e [a,b]) is minimal, i.e. dim Nt{S^'''^'^x) = 1. 

(2) The trajectory segment S^"''^^x containing exactly one singularity (a so called 
"simple singularity" , see above) is said to be sufficient if and only if both 
branches of this trajectory segment are sufficient. 

Definition 2.8. The phase point a; e M with at most one (simple) singularity is 

said to be sufficient if and only if its whole trajectory S^~°°'°°'^x is sufficient, which 
means, by definition, that some of its bounded segments S^°''''^x are sufficient. 
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Note. In this paper the phrase "trajectory (segment) with at most one singularity" 
always means that the sole singularity of the trajectory (segment), if exists, is simple. 

In the case of an orbit S^~°°'°°^x with at most one singularity, sufficiency means 
that both branches of S^~°°'°°^x are sufficient. 

No accumulation (of collisions) in finite time. By the results of Vaserstein 
[V(1979)], Galperin [G(1981)] and Burago-Ferleger-Kononenko [B-F-K(1998)], in 

any semi-dispersive billiard flow there can only be finitely many collisions in finite 
time intervals, see Theorem 1 in [B-F-K(1998)]. Thus, the dynamics is well defined 
as long as the trajectory does not hit more than one boundary components at the 
same time. 

Slim sets. We are going to summarize the basic properties of codimension-two 
subsets A of a connected, smooth manifold M with a possible boundary and corners. 
Since these subsets A are just those negligible in our dynamical discussions, we shall 
call them slim. As to a broader exposition of the issues, see [E(1978)] or §2 of [K- 
S-Sz(1991)]. 

Note that the dimension dim A of a separable metric space A is one of the three 
classical notions of topological dimension: the covering (Cech-Lebesgue), the small 
inductive (Menger-Urysohn), or the large inductive (Brouwer-Cech) dimension. As 
it is known from general topology, all of them are the same for separable metric 
spaces, see [E(1978)]. 

Definition 2.9. A subset ^ of M is called slim if and only if A can be covered by 
a countable family of codimension-two (i. e. at least two) closed sets of |U-measure 
zero, where n is any smooth measure on M. (Cf. Definition 2.12 of [K-S-Sz(1991)].) 

Property 2.10. The collection of all slim subsets of M is a cr-ideal, that is, count- 
able unions of slim sets and arbitrary subsets of slim sets are also slim. 

Proposition 2.11. (Locality). A subset A G M is slim if and only if for every 
X E A there exists an open neighborhood U of x in M such that U DA is slim. (Cf. 
Lemma 2.14 of [K-S-Sz(1991)].) 

Property 2.12. A closed subset ^ C M is slim if and only if //(A) = and 
dim^ < dimM - 2. 

Property 2.13. (Integrability) . If ^ C Mi x M2 is a closed subset of the product 
of two smooth, connected manifolds with possible boundaries and corners, and for 
every x e Mi the set 

A^ = {ye M2: {x,y) e A} 

is slim in M2, then A is slim in Mi x M2. 

The following propositions characterize the codimension-one and codimension- 
two sets. 
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Proposition 2.14. For any closed subset S G M the following three conditfons 
are equivalent: 

(i) dim 5 < dimM - 2; 

(ii) int5' = and for every open connected set G C M the difference set G\S 
is also connected; 

(iii) intS' = and for every point x e M and for any open neighborhood V of x 
in M there exists a smaller open neighborhood W CV of the point x such 
that for every pair of points y, z & W \ S there is a continuous curve 7 in 
the set V\S connecting the points y and z. 

(See Theorem 1.8.13 and Problem 1.8.E of [E(1978)].) 

Proposition 2.15. For any subset S C M the condition dim 5' < dimM — 1 is 
equivalent to int^ = 0. (See Theorem 1.8.10 of [E(1978)].) 

We recall an elementary, but important lemma (Lemma 4.15 of [K-S-Sz(1991)]). 
Let A2 be the set of phase points a; G M \ 9M such that the trajectory S^~°°'°°^x 
has more than one singularities (or, its only singularity is not simple). 

Proposition 2.16. The set A2 is a countable union of codimension-two smooth 
submanifolds of M and, being such, is slim. 

The next lemma establishes the most important property of slim sets which 
gives us the fundamental geometric tool to connect the open ergodic components 
of billiard flows. 

Proposition 2.17. If M is connected, then the complement M\A of a slim Fcr set 
A C M is an arc- wise connected {Gg) set of full measure. (See Property 3 of §4.1 in 
[K-S-Sz(1989)]. The sets are, by definition, the countable unions of closed sets, 
while the G5 sets are the countable intersections of open sets.) 

The subsets M*^ and M#. Denote by M# the set of all phase points x G M for 
which the trajectory of x encounters infinitely many non-tangential collisions in both 
time directions. The trajectories of the points a; G M \ are lines: the motion is 
linear and uniform, see the appendix of [Sz(1994)]. It is proven in lemmas A. 2.1 and 
A. 2. 2 of [Sz(1994)] that the closed set M \ is a finite union of hyperplanes. It 
is also proven in [Sz(1994)] that, locally, the two sides of a hyper-planar component 
of M \ M'^ can be connected by a positively measured beam of trajectories, hence, 
from the point of view of ergodicity, in this paper it is enough to show that the 
connected components of entirely belong to one ergodic component. This is 
what we are going to do in this paper. 

Denote by M*^ the set of all phase points x G M# the trajectory of which does not 
hit any singularity, and use the notation for the set of all phase points x G M'^ 
whose orbit contains exactly one, simple singularity. According to Proposition 
2.16, the set M# \ (M" U M^) is a countable union of smooth, codimension-two 
(> 2) submanifolds of M, and, therefore, this set may be discarded in our study of 
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ergodicity, please see also the properties of slim sets above. Thus, we will restrict 
our attention to the phase points x e M° U M^. 

The "Chernov-Sinai Ansatz". An essential precondition for the Theorem on 
Local Ergodicity by Chernov and Sinai [S-Ch(1987)] is the so called "Chernov- 
Sinai Ansatz" which we are going to formulate below. Denote by STZ^ C the 
set of all phase points xq = (90,^0) ^ dNl corresponding to singular reflections 
(a tangential or a double collision at time zero) supplied with the post-collision 
(outgoing) velocity vq. It is well known that STZ'^ is a compact cell complex with 
dimension 2d — 3 = dimM — 2. It is also known (see Lemma 4.1 in [K-S-Sz(1990)-I], 
in conjunction with Proposition 2.16 above) that for i^i-almost every phase point 
Xq G STZ'^ the forward orbit S^'^'^^^xq does not hit any further singularity. (Here 
1^1 is the Riemannian volume of STZ'^ induced by the restriction of the natural 
Riemannian metric of M.) The Chernov-Sinai Ansatz postulates that for i/i-almost 
every a^o £ STZ'^ the forward orbit S^^'°°''xo is suflBcient (geometrically hyperbolic). 

The Theorem on Local Ergodicity. The Theorem on Local Ergodicity for 
semi-dispersive billiards (Theorem 5 of [S-Ch(1987)]) claims the following: Let 
(M, {jS^lfgiK, /u) be a semi-dispersive billiard flow with the property that the smooth 
components of the boundary dQ, of the configuration space are algebraic hyper- 
surfaces. (The cylindric billiards automatically fulfill this algebraicity condition.) 
Assume - further - that the Chernov-Sinai Ansatz holds true, and a phase point 
Xq e (M° U M^) \ aM is sufficient. 

Then some open neighborhood [/q C M of xq belongs to a single ergodic compo- 
nent of the flow (M, {5'*}teK5 A*)- (Modulo the zero sets, of course.) 

53. Expansion and Contraction Rate Estimates 



We would like to get a useful lower estimate for the expansion of a tangent vec- 
tor {Sqo,dvo) G M with positive infinitesimal Lyapunov function Q{5qo,Svo) = 
{5qo,5vo). The expression {5qo,5vo) is the scalar product in R.^ defined via the 
mass (or kinetic energy) metric, see §2. It is also called the infinitesimal Lyapunov 
function associated with the tangent vector (dqa^dvo), see [K-B(1994)], or part A. 4 
of the Appendix in [Ch(1994)], or §7 of [Sim(2003)]. For a detailed exposition of the 
relationship between the quadratic form Q{.), the relevant symplectic geometry of 
the Hamiltonian system and the dynamics, please also see [L-W(1995)]. 

Note. The original idea of using infinitesimal Lyapunov exponents to measure 
the expansion rate of codimension-one submanifolds in the phase space of semi- 
dispersive billiards came from N. Chernov back in the late '80s. These ideas have 
been explored in detail and further developed by him and myself in recent personal 
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communications, so that we obtained at least linear (but uniform!) expansion rates 
for such submanifolds with negative infinitesimal Lyapunov exponents for their 
normal vector. These results will be presented in our forthcoming joint paper [Ch- 
S(2006)], and they will be used later in this article in the proof of Lemma 6.2. Also, 
closely related to the above said, the following ideas (to estimate the expansion rates 
of tangent vectors from below) are derived from the thoughts mentioned above. 

Denote by {Sqt, 6vt) — {DS^){Sqo, Svq) the image of the tangent vector {Sqo, Svq) 
under the linearization DS* of the map S*, t >0. (We assume that the base phase 
point xq — for which (Sqo, Svq) e — has a non-singular forward orbit.) The 

time-evolution {Sqtj^jdvt^) ^ [Sqt^^Svt^] (0 < ti < t2) on a collision free segment 
5'[ti,t2]2>Q is described by the equations 

(3 1) ^'"t.^Svt,, 

Correspondingly, the change Q{Sqti, Svt^) i— > Q{5qt^,5vt2) in the infinitesimal Lya- 
punov function Q{.) on the collision free orbit segment S^^'^'^'^^xq is 



(3.2) Q{dqt,,dvt,) = Q{dqt,,Svt,) + {t2 - h)\\dvt,\\^ 

thus Q{ .) steadily increases between collisions. 

The passage {5q^,Svl^) h- >• {6q^,5v^) through a reflection (i. e. when xt = 
S^xo e dM) is given by Lemma 2 of [Sin(1979)] or formula (2) in §3 of [S-Ch(1987)]: 

5q^ = RSqr, 

(3.3) * ' 

5vt = R5vt + 2 cos 4>RV*KVSqt , 

where the operator R : TCI ~^ is the orthogonal reflection (with respect to 
the mass metric) across the tangent hyperplane Tg^dQ, of the boundary dQ, at the 
configuration component qt of Xt = {qt,vf), V : (v^)-^ Tq^dQ is the t;^~-parallel 
projection of the orthocomplement hyperplane (ft )"*" onto Tq^dQ, V* : Tq^dQ —>■ 
(v^)-^ is the adjoint of V (i. e. the z>'(gt)-parallel projection of Tq^dQ onto (ft")"*", 
where I'iqt) is the inner normal vector of dQ at qt G 5Q), K : Tq^dQ — > Tq^dQ is 
the second fundamental form of the boundary dQ, at qt (with respect to the field 
^{q) of inner unit normal vectors of dQ) and, finally, coscp — {y{qt), v^) > is 
the cosine of the angle (0 < < 7r/2) subtended by and z/(gt). Regarding 
formulas (3.3), please see the last displayed formula in §1 of [S-Ch(1987)] or (i)-(ii) 
in Proposition 2.3 of [K-S-Sz(1990)-I]. The instanteneous change in the infinitesimal 
Lyapunov function Q{Sqt, Svt) caused by the refiection at time t > is easily derived 
from (3.3): 
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Q{Sq+,Sv+) = Q{Sqi,Sv^) + 2coscl>{VSqi,KVSqi) 

(3.4) 

> (3(5?^ ,Svt ). 

In the last inequality we used the fact that the operator K is positive semi-definite, 
i. e. the billiard is semi-dispersive. 

We are primarily interested in getting useful lower estimates for the expansion 
rate ||5q't||/||5q'o||- The needed result is 

Proposition 3.5. Use all the notations above, and assume that 

{6qo, Svo)/\\6qo\\^ > Co > 0. 
We claim that | I^Qt 1 1/| |5qo| I > 1 + cot for all t>0. 

Proof. Clearly, the function \ \6qt\ \ of t is continuous for alH > and continuously 
differentiable between collisions. According to (3.1), 4:Sqt — Svt-, so 



(3.6) -\\5qt\\^ = 2{5qt,5vt). 

Observe that not only the positive valued function Q{5qt, 6vt) = {Sqt,Svt) is 
nondecreasing in t by (3.2) and (3.4), but the quantity {Sqt, Svt) /\\dqt\\ is nonde- 
creasing in t, as well. The reason is that {Sqt, Svt) /\\dqt\\ = \\dvt\\cosat {at being 
the acute angle subtended by 6qt and dvt), and between collisions the quantity \ \Svt\\ 
is unchanged, while the acute angle at decreases, according to the time-evolution 
equations (3.1). Finally, we should keep in mind that at a collision the norm \ \Sqt\\ 
does not change, while {Sqt,Svt) cannot decrease, see (3.4). Thus we obtain the 
inequalities 



{Sqt,dvt)/\\dqt\\ > {dqo, dvo) /\\dqo\\ > Co\\dqo\\, 

so 

f^\\Sqt\f = 2\\6qt\\f^\\6qt\\ = 2{6qt,6vt) > 2co\\6qo\\ ■ \\6qt\\ 

by (3.6). This means that ^11(^9*11 > coH^g-oH, so \\6qt\\ > \\6qo\\{l + cot) , proving 
the proposition. □ 

Next we need an efi^ective lower estimation cq for the curvature {dqo, dvo) /\\Sqo\\'^ 
of the trajectory bundle: 
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Lemma 3.7. Assume that the perturbation {Sqq , Svq) e %^M. (as in Proposition 
3.5) is being performed at time zero right before a coUision, say, ctq = (1, 2) taking 
place at that time. Select the tangent vector {Sqq , Svq) in such a specific way that 
Svq = 0, Sqq = {m2W, —miw, 0, 0, . . . , 0) with a nonzero vector w G M'^, {w, — 
V2) = 0. This scalar product equation is exactly the condition that guarantees that 
Sqq be orthogonal to the velocity component v~ — {v^ ^ V2 , . . . , v^) of xq — {q, v~). 
The last, though crucial requirement is that w should be selected from the two- 
dimensional plane spanned by f ^ — f ^ and qi — q2 (with | l^i — q'2 1 1 = 2r) in M'^. The 
purpose of this condition is to avoid the unwanted phenomenon of "astigmatism" in 
our billiard system, discovered first by Bunimovich and Rehacek in [B-R(1997)] and 
[B-R(1998)]. Later on the phenomenon of astigmatism gathered further prominence 
in the paper [B-Ch-Sz-T(2002)] as the main driving mechanism behind the wild 
non-diff'erentiability of the singularity manifolds (at their boundaries) in hard ball 
systems in dimensions bigger than 2. We claim that 



.0 o^ ^ ^•^0 ) 11^1 -'^2! I 11^1 --^2|| 

||()go|r r cos 00 r 

for the post-collision tangent vector {5qQ,5vQ), where 00 is the acute angle sub- 
tended by Vi — V2 and the outer normal vector of the sphere {y e ||y|| = 2r} 
at the point y = qi — q2- Note that in (3.8) there is no need to use + or — in ||5go|P 
or \\vi - V2\\, for \\Sq^\\ = \\Sq^\\, \\v~ - v^W = \\vf -v^\\. 

Proof. The proof of the equation in (3.8) is a simple, elementary geometric argu- 
ment in the plane spanned by — V2 and qi — q2, so we omit it. We only note 
that the outgoing relative velocity f ^ — is obtained from the pre-collision relative 
velocity f f — by reflecting the latter one across the tangent hyperplane of the 
sphere {y e M^^l ||j/|| = 2r} at the point y = qi — q2- D 

The previous lemma shows that, in order to get useful lower estimations for the 

"curvature" {Sq,Sv) /WSqW^ of the trajectory bundle, it is necessary (and sufficient) 
to find collisions a = on the orbit of a given point xq G M with a "relatively 
big" value of — Finding such collisions will be based upon the following 
result: 

Proposition 3.9. Consider orbit segments S^^'^^xq of A'^-ball systems with masses 
mi, 1712, ■ ■ ■ ,mN in T'^ (or in M^) with collision sequences E = (ui, a"2, . . . , cr„) cor- 
responding to connected collision graphs. (Now the kinetic energy is not necessarily 
normalized, and the total momentum Y2iLi '^i'^i may be different from zero.) We 
claim that there exists a positive- valued function /(a; mi, m2, . . . , m^r) (a > 0, / is 
independent of the orbit segments S^^^'^^Xq) with the following two properties: 

(1) If ||fi(tO — fj(^z)ll — collisions ct; = (z,j) (1 < / < n., is the time 

of Gi) of some trajectory segment S^^^'^^xq with a symbolic collision sequence E = 
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(o"! , (72, . . . , corresponding to a connected collision graph, then the norm | |f j/ (t) — 
Vj'{t)\ \ of any relative velocity at any time t e R is at most /(a; mi, . . . jTUn)', 
(2) lima^o /(a; mi, ... , uin) = for any (mi, ... , uin)- 

Proof. We begin with 

Lemma 3.10. Consider an AT-ball system with masses mi, . . . ,mjv (an (mi, . . . 
,mjv)-system, for short) in (or in R'^). Assume that the inequalities ||vi(0) — 
Vj(0)|| < a hold true {I < i < j < N) for all relative velocities at time zero. We 
claim that 



(3.11) \\vi{t) -Vj{t)\\ <2a 

for any pair (i, j) and any time t G R, where M = X^^L^ m^, 

m = min {rrii \ 1 < i < N} . 

Note. The estimate (3.11) is far from the optimal one, however, it will be sufficient 
for our purposes. 

Proof. The assumed inequalities directly imply that ||fi(0)|| < a (1 < z < N) 
for the velocities ^^(0) measured at time zero in the baricentric reference system. 
Therefore, for the total kinetic energy Eq (measured in the baricentric system) we 
get the upper estimation Eq < ^Ma^, and this inequality remains true at any time 
t. This means that all the inequalities ||v^(t)|P < ^cb"^ hold true for the baricentric 
velocities vl{t) at any time t, so 

\\vl{t) - v'j{t)\\ < aVM (mr'/' + m/'/') < 2a 
thus the inequalities 

\\vi{t) -Vj{t)\ \ < 2a 

hold true, as well. □ 

Proof of the proposition by induction on the number N. 

For = 1 we can take /(a; mi) = 0, and for A = 2 the function /(a; mi, 7712) — a 
is obviously a good choice for /. Let A > 3, and assume that the orbit segment 
S^^''^^xo of an (mi, . . . , mjv)-system fulfills the conditions of the proposition. Let 
CTfc = be the collision in the symbolic sequence E„ = (cti, . . . , fT^j,) of S^'^^'^^xq 

with the property that the collision graph of = (a"i, . . . , 0"^) is connected, while 
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the collision graph of — (cti, . . . , crfe_i) is still disconnected. Denote the two 

connected components (as vertex sets) of by Ci and C2, so that i e Ci, 

j e C2, C1UC2 = {1,2,..., N}, and Ci n C2 = 0- By the induction hypothesis 
and the condition of the proposition, the norm of any relative velocity Vi'{tk — 
0) — Vj'{tk — 0) right before the collision Ufe (taking place at time tk) is at most 
a + f{a; Ci) + f{a] C2), where C/ stands for the collection of the masses of all 
particles in the component Ci, I = 1, 2. Let g{a; mi, . . . , mjv) be the maximum of 
all possible sums 

a + f{a; Di) + f{a;D2), 

taken for all two-class partitions (Di, D2) of the vertex set {1,2,..., N}. According 
to the previous lemma, the function 

IW 

f{a; mi, ... , uin) := 2W — g{a; mi, ... , mjv) 
fulfills both requirements (1) and (2) of the proposition. □ 

Corollary 3.12. Consider the original (mi, . . . , mAr)-system with the standard 

normalizations X^ili "^i'^^i ~ ^ '^iLi = ^- We claim that there exists a 

threshold G = G{mi, . . . jUIn) > (depending only on N, mi, . . . jUIn) with the 
following property: 

In any orbit segment S^^''^^xo of the (mi, . . . , mAr)-system with the standard 
normalizations and with a connected collision graph, one can always find a collision 
cr = taking place at time t, so that ||fi(t) — 't^j(^)|| > G(mi, . . . jUIn)- 

Proof. Indeed, we choose G = G{mi, . . . ,mN) > so small that f{G; mi, . . . , mjv) 
< M~^/'^. Assume, contrary to 3.12, that the norm of any relative velocity Vi — Vj 
of any collision of S^^''^'^xo is less than the above selected value of G. By the 
proposition, we have the inequalities ||'i'i(0) — ^^(O)!! < /(G; mi, . . . , m^v) at time 
zero. The normalization ^^iTOifi(O) = 0, with a simple convexity argument, 
implies that ||vi(0)|| < f{G; mi, . . . ,mjv) for alH, 1 < z < N, so the total kinetic 

1 2 1 

energy is at most [f{G; mi, . . . , mjv)] < 5, a contradiction. □ 

Corollary 3.13. For any phase point xq with a non-singular backward trajec- 
tory S^~°^'^^xo and with infinitely many consecutive, connected collision graphs 

on S^~°°'^^Xo, and for any munber L > one can always find a time —t<0 
and a tangent vector {dqo,Svo) G T^ jM {x-t = S~^xo) with {Sqo,Svo) > and 
ll^^tll/ll^^oll > L, where {5qt,5vt) ^ DS\5qo,5vo). 

Proof. Indeed, select a number t > so big that l-|-^G(mi, . . . jUIn) > L and — t is 
the time of a collision (on the orbit of xq) with the relative velocity v~{—t) —v~{—t), 
for which \\v^{—t) — Vj~{—t)\\ > G(mi, . . . , mjv)- By Lemma 3.7 we can choose a 
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tangent vector [dq^ , 0) right before the coUisfon at time —t in such a way that the 
fower estimation 

{dq+, dv+) , 1 s 

holds true for the "curvature" {Sqq , ^'^o) /W^^oW^ associated with the post-colhsion 
tangent vector {SqQ, Svq). According to Proposition 3.5 we have then the fower 
estimation 

||M>i + !G(mi,...,mAr)>L 
for the (Jq'-expansion rate between {Sqq , 0) and {5qt, 5vt) = DS*{5qQ, 0). □ 



§4. The Exceptional J-Manifolds 
(The asymptotic measure estimates) 



First of all, we define the fundamental object for the proof of our theorem. 

Definition 4.1. A smooth submanifold J C intM of the interior of the phase space 
M is called an exceptional J -manifold (or simply an exceptional manifold) with a 
negative Lyapunov function Q if 

(1) dimJ = 2d - 2 (= dimM - 1); 

(2) the pair of manifolds (J, dJ) is diffeomorphic to the standard pair 

where B^^~'^ is the closed unit ball of E^'^"^; 

(3) J is locally flow-invariant, i. e. \/x E J 3a{x), b{x), a{x) < < b{x), such 
that S*x e J for aU t with a{x)<t< b{x), and 5"(^)a; e dJ, S^'^'^^x e dJ; 

(4) the manifold J has some thin, open, tubular neighborhood Uq in intM, and 
there exists a number T > such that 

(i) j n dM. = 0, and all orbit segments S't^'^^x {x e Uq) are non-singular, 
hence they share the same symbolic collision sequence E; 

(ii) Vx e JJq the orbit segment S^^''^^x is sufficient if and only if a; ^ J; 

(5) Vx e J we have Q{n{x)) := {z{x), w{x)) < — ci < for a unit normal vector 
field n{x) = {z{x), w{x)) of J with a fixed constant ci > 0; 
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(6) the set W of phase points x E J never again returning to J (After first leaving 
it, of course. Keep in mind that J is locally flow-invariant!) has relative measure 

greater than 1 — 10~^ in J, i. e. > 1 — 10~^, where /ii is the hypersurface 

^i(J) 

measure of the smooth manifold J. 



We begin with an important proposition on the structure of forward orbits 

Proposition 4.2. For //i-almost every x E J the forward orbit S^^'°°''x is non- 
singular. 

Proof. According to Proposition 7.12 of [Sim(2003)], the set 



Jn 



U i^^-) 



.t>0 



of forward singular points a; G J is a countable union of smooth, proper submanifolds 
of J, hence it has //i-measure zero. □ 

In the future we will need 

Lemma 4.3. The concave, local orthogonal manifolds S(y) passing through points 
y E J are uniformly transversal to J. 

Note. A local orthogonal manifold E C intM is obtained from a codimension-one, 
smooth sub manifold Ei of intQ by supplying Ei with a selected field of unit normal 
vectors as velocities. E is said to be concave if the second fundamental form of Ei 
(with respect to the selected orientation) is negative semi-definite at every point of 
El. Similarly, the convexity of S requires positive semi-definiteness here, see also 
§2 of [K-S-Sz(1990)-I]. 

Proof. We will only prove the transversality. It will be clear from the uniformity 
of the estimations used in the proof that the claimed transversalities are actually 
uniform across J. 

Assume, to the contrary of the transversality, that a concave, local orthogonal 
manifold E(y) is tangent to J at some y & J. Let {Sq, B5q) be any vector of TyM. 
tangent to E(y) at y. Here S < is the second fundamental form of the projection 
g(E(y)) = Ei(y) of E(y) at the point q = q{y). The assumed tangency means that 
{dq, z) + (Bbq^ w) — 0, where n{y) = {z{y), w{y)) = {z, w) is the unit normal vector 
of J at y. We get that {6q, z + Bw) — for any vector Sq e v{y)-^. We note that the 
components z and w of n are necessarily orthogonal to the velocity v{y), because the 
manifold J is locally flow-invariant. The last equation means that z = —Bw, thus 
Q{n{y)) = {z, w) = {—Bw, w) > 0, contradicting to the assumption Q{n{y)) < —c\ 
of (5) in 4.1. This flnishes the proof of the lemma. □ 
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In order to formulate the main result of this section, we need to define two 
important subsets of J. 

Definition 4.4. Let 

A=^xeJ\ S^^'°°^x is nonsingular and dimjVo (^S^^'°°^x^ = l} , 

B = ^xeJ\ S^'^'°°^x is nonsingular and dimjVo (^S^^'°°^x^ > l} . 

The two Borel subsets A and 5 of J are disjoint and, according to Proposition 
4.2 above, their union AU B has full //i-measure in J. 
The anticipated main result of this section is 

Main Lemma 4.5. Use all of the above definitions and notations. We claim that 

Proof. The proof will be a proof by contradiction, and it will be subdivided into 
several lemmas. Thus, from now on, we assume that A = (/). 

First, select and fix a non-periodic point (a "base point") xq G B. The following 
step is to use Corollary 3.13 in such a way that the role of xq in 3.13 is now 
played by the time-inverted version —xq = (go, —vq) of our fixed base point xq e J. 
Thus, for a large constant Lq » 1 (to be specified later) select a big enough time 
C3 >> 1 (playing the role of t in 3.13) and a tangent vector {Sqo, —Svq) e Tl^j^^M 
(xcg = S'^^xq) with {Sqo, —Svq) > and 



(4.6) ^ > L„, 

ll^^oll 

where 



(4.7) (Sqc,, Svc,) = {DS-'^') {Sqo, Svo). 

The normalized tangent vector 



(4.8) {Sqo, Svo) ■■= (I|5gc3ll' + ll^^^call') '^^-{Sqcs. Sv,,) G %,M 

will play a crucial role in the proof. 

By slightly perturbing the tangent vector {Sqo, —Svo), we can always achieve that 



the unit tangent vector {Sqo, Svq) of (4.8) 
(4 9) ) is transversal to J and, corresponding to (4.6), 
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where {Sqc^, dvcs) = {DS''^) {dqo, Svq). 

We choose the orientation of the unit normal field n{x) {x e J) of J in such a 
way that {n{xo), {dqo, Svq)) < 0, and define the one-sided tubular neighborhood Us 
of radius 5 > as the set of all phase points 7x(s), where x G J, < s < 5. Here 
"Yxi-) is the geodesic line passing through x (at time zero) with the initial velocity 
n{x), X e J. The radius (thickness) (5 > here is a variable, which will eventually 
tend to zero. We are interested in getting useful asymptotic estimates for certain 
subsets of Us, as 5 — > 0. 

Our main working domain will be the set 

Dq^ & Us^\j\y S~* {Sn~) , 3 a sequence 
(4.10) *>o 

tn y oo such that S^'^y eUso\J, n = 1, 2, . . . |, 

a set of full |U-measure in Us^- We will use the shorthand notation Uq = Us^ for a 
fixed, small value 5q of 5. For any y e M we use the traditional notations 



riy) = min {t > 0| S^y e dM] , 

(4.11) ' ^ 

T{y) = S^y^y 

for the first hitting of the collision space 9M. The first return map (Poincare 
section, collision map) T : dM. — > dM. (the restriction of the above T to dM) is 
known to preserve the finite measure f that can be obtained from the Liouville 
measure ^ by projecting the latter one onto dM along the fiow. Following 4. of 
[K-S-Sz(1990)-II], for any point y G intM (with r(y) < oo, T{—y) < oo, where 
~y — {Q: ~'>j) ^oi y = {q, v)) we denote by Ztub{y) the supremum of all radii p > 
of tubular neighborhoods Vp of the projected segment 

q{{S'y\ -T(-y)<t<T(y)}) CQ 
for which even the closure of the set 



{(g, v{y))&M\qeVp] 

does not intersect the set STZ of singular reflections. 

We remind the reader that both Lemma 2 of [S-Ch(1987)] and Lemma 4.10 of 
[K-S-Sz(1990)-I] use this tubular distance function Ztub{-) (despite the notation 
2;( . ) in those papers), see the important note 4. on [K-S-Sz(1990)-II]. 

Following the fundamental construction of local stable invariant manifolds [S- 
Ch(1987)] (see also §5 of [K-S-Sz(1990)-I]), for any y e Dq we define the concave, 
local orthogonal manifolds 
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^l{y) = SCy, {{{q, vivt)) e M| g - g(y,) ± v{yt)} \ {S, U , 
= SCy [S-'^l{y)] , 

where iSi := {x G M| Tx G 57^~} (the set of phase points on singularities of order 
1), S-i := {a; G M| — x E Si^ (the set of phase points on singularities of order 
— 1), ut — S^y, and SCy{ . ) stands for taking the smooth component of the given 
set that contains the point y. The local, stable invariant manifold 7^*^(y) of y is 
known to be a superset of the C^-limiting manifold limt^oo ^oiu)- 

On all these local orthogonal manifolds, appearing in the proof, we will always 
use the so called 5g-metric to measure distances. The length of a smooth curve 
with respect to this metric is the integral of | |5g| | along the curve. The proof of the 
Theorem on Local Ergodicity [S-Ch(1987)] shows that the (Jg-metric is the relevant 
notion of distance on the local orthogonal manifolds E, also being in good harmony 
with the tubular distance function Ztub{ ■ ) defined earlier. 

On any manifold Sgd/) fl C/q [v ^ Dq) we define the smooth field Xy^t{y') {y' G 
Sfld/) n Uq) of unit tangent vectors of So(2/) ^ as follows: 



(4.13) Xy,iy')= n.,t„'((5go,5^)o)) 



where Iiy,t,y' denotes the orthogonal projection of M.^ © R'' onto the tangent space 
of So(y) at the point y' G So(y) H Uq. Recall that {Sqo, Svq) is the unit tangent 
vector of M at the base point xq from (4.8) and (4.9). We also remind the reader 
that {Sqo, Svq) points toward the side of J opposite to the side where the one-sided 
neighborhoods U5 reside. 

Note 4.14. By the construction of {6qo,Svo) in (4.6)-(4.8), if the threshold C3 is 
big enough, then the vector {dqo, Svq) is close to the tangent space T'y^{xo) of the 
local stable manifold of a^o- On the other hand, for large enough t the tangent space 
of Sq^/) n Uq at y' makes a small angle with T7*(a;o). All the necessary upper 
estimations for the mentioned angles follow from the well known result stating that 
the difference (in norm) between the second fundamental forms of the 5'*-images 
{t > 0) of two local, convex orthogonal manifolds is at most see, for instance, 
inequality (4) in [Ch(1982)]. These facts imply that the vector in the numerator of 
(4.13) is actually very close to {5qo, Svq), in particular its magnitude is almost one. 

For any y e Dq let tk = tk{y) (0 < ti < t2 < • • • ) be the time of the A;-th collision 
CTfe on the forward orbit S^^'°°^y of y. Assume that the time t in the construction 

of So(?/) and Xy^t is between (Tk-i and CTfc, i. e. tk-i{y) < t < tk{y)- We define the 
smooth curve py^t = Py,t{.s) (with the arc length parametrization s, < s < h{y,t)) 
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as the maximal integral curve of the vector field Xy t emanating from y and not 
intersecting any forward singularity of order < k, i. e. 

Py,t(0) = y, 

Py,t{ ■ ) does not intersect any singularity of order < fc, 
Py^t is maximal among all curves with the above properties. 

We remind the reader that a phase point x lies on a singularity of order (fc e M) 
if and only if the fc-th collision on the forward orbit S'^'^'°°^x is a singular one. It 
is also worth noting here that, as it immediately follows from (4.15), the curve py^t 
can only terminate at a boundary point of the manifold TiQ^y) H Uq. 

Note 4.16. From now on, we will use the notations So(y), Xy^^, and Py^^ foi" 

^o{y)^ ^y,ti, and py,t., respectively, where tl = tl{y) = ^ {tk-i{y) +tk{y)). 

Due to these circumstances, the curves Py,ti — Py,k can now terminate at a point 
z such that z is not on any singularity of order at most fc and S^^z is a boundary 
point of Th (y) , so that at the point S^^ z the manifold E*J (y) touches the boundary 

k k 

of the phase space in a nonsingular way. This means that, when we continuously 
move the points Py,k{s) by varying the parameter s between and h{y, fc), either 
the time {Py,k{s)) or the time t^-i (Py.ki-^)) becomes equal to t^. = tl.{y) when the 
parameter value s reaches its maximum value h{y, fc). The length of the curve Py^k is 
at most So, and an elementary geometric argument shows that the time of collision 
tk{py,k{s)) (or tk-i{py,k{s))) can only change by at most the amount of c*-\/5o, as 
s varies between and h{y,k). (Here c* is an absolute constant.) Thus, we get 
that the unpleasant situation mentioned above can only occur when the difference 
tk{y) — tk-i{y) is at most c*y/So. These collisions have to be excluded from the 
limiting process of Lemma 4.28 below. Still, everything works by the main result of 
[B-F-K(1998)], which guarantees that the indices fc of the collisions (on the forward 
orbit S^^'°°^y) with tk{y) — tk-i{y) > c*\/5q have a positive density amongst the 
natural numbers. 

As far as the terminal point Py^k{h{y, k)) of Py^k is concerned, there are exactly 
three, mutually exclusive possiblities for this point: 

(A) py^kih{y,k)) e J and this terminal point does not belong to any forward 
singularity of order < fc, 

(B) py^kih{y, fc)) lies on a forward singularity of order < fc, 

(C) the terminal point py^k{h{y, fc)) does not lie on any singularity of order < fc 
but lies on the part of the boundary dUo of Uq difi^erent from J. 
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Note 4.17. Under the canonical identification Uq = J x [0, do) of Uq via the 
geodesic fines perpendicular to J, the above mentioned part of dUo (the "side" of 
Uq) corresponds to dJx [0, 6o). Therefore, the set of points with property (C) inside 
a layer Us {S < do) will have |U-measure small ordo of d (actually, of order 5^), thus 
this set will be negligible in our asymptotic measure estimations, as 5 ^ 0. That is 
why in the future we will not be dealing with any phase point with property (C). 

Should (B) occur for some value of k {k > 2), the minimum of all such integers 
k will be denoted hj k = k{y). The exact order of the forward singularity on which 
the terminal point ^ {h{y, k)) lies is denoted by ki — ki{y). If (B) does not occur 

for any value of k, then we take k{y) — ki{y) = oo. 

We can assume that the manifold J and its one-sided tubular neighborhood 
Uq = are already so small that for any y E Uq no singularity of S^^'"°^y can 
take place at the first collision, so the indices k and ki above are automatically at 
least 2. For our purposes the important index will be ki — ki{y) for phase points 
y e i^o- 

Note 4.18. Refinement of the construction. Instead of selecting a single 
contracting unit vector {6%, Svq) in (4.8), we should do the following: Choose a 
compact set Kq C B with the following properties: 

(i) ^ > 1 - lo-^ 

(ii) every point x e Kq has a non-singular forward orbit. 

Now the running point x G Kq will play the role of xq in the construction of the 

contracting unit tangent vector u{x) := {5%, Svq) G T^M on the left-hand-side of 
(4.8). For every x G Kq there is a small, open ball neighborhood B{x) of x and a 
big threshold C3(x) >> 1 such that (4.9) holds true for u{x) and C3 = 03(0;) for all 

X G Kq. 

By the continuity of the contraction/expansion factor, one can also achieve that 
the contraction estimation Lq^ of (4.9) holds true not only for u{x), but also for 
any projected copy of it appearing in (4.13), provided that y' G B{x), i. e. y' is 
close enough to x. 

Now select a finite subcovering IJ^^i B{xi) of Kq, and replace J by Ji = J fl 
[J^^i B{xi), Us by U'g = Us n Ur=i -^(^») (^^^ ^ — ^0) finally, choose the 

threshold C3 to be the maximum of all thresholds cs(xi) for z = 1, 2, . . . , n. In this 
way the assertion of Corollary 4.20 will be indeed true. 

We note that the new exceptional manifold Ji is no longer so nicely "roimd 
shaped" as J, but it is still pretty well behaved, being a domain in J with a piecewise 
smooth boundary. 

The reason why we cannot switch completely to a round and much smaller man- 
ifold B{x) n J is that the measure fii{J) should be kept bounded from below after 
having fixed Lq, see 4. in the Appendix. 
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In addition, it should be noted that, when constructing the vector field in (4.13) 
and the curves py^ti an appropriate directing vector u{xi) needs to be chosen for 
(4.13). To be definite and not arbitrary, a convenient choice is the first index 
i e {1, 2, . . . , n} for which y G B{xi). In that way the whole curve py^t will stay 
in the slightly enlarged ball B'{xi) with double the radius of B{xi), and one can 
organize things so that the required contraction estimates of (4.9) be still valid even 
in these enlarged balls. 

In the future, a bit sloppily, Ji will be denoted by J, and U'g by Us- 

Note 4.19. When defining the returns of a forward orbit to Us, we used to say 
that "before every new return the orbit must first leave the set Us"- Since the 
newly obtained J is no longer round shaped as it used to be, the above phrase is 
not satisfactory any longer. Instead, one should say that the orbit leaves even the 
K-neighborhood of Us, where k, is two times the diameter of the original J. This 
guarantees that not only the new Us, but also the original Us will be left by the 
orbit, so we indeed are dealing with a genuine return. This note also applies to two 
more slight shrinkings of J that will take place later in the proof. 

As an immediate corollary of (4.9) and the above note, we get 

Corollary 4.20. For the given sets J, Uq, and the large constant Lq we can select 
the threshold C3 > large enough so that for any point y G Dq any time t with 
C3 < t < 5q'-expansion rate of 5"* between the curves Py^y^ and 

{Pyjii^y)j is less than Lq \ i. e. for any tangent vector {5qQ, 5vo) of Py^^^y-^ we 
have 



where {5qt, Svt) = {DS*){6qo, 6vo)- 

A further immediate consequence of the previous result is 

Corollary 4.21. For any y e Dq with k{y) < 00 and t^_^^y^_i{y) > C3, and for any 
t with %^(y)_i(y) < t < %^y){y), we have 

(4.22) ztub {S'y) < Lq\ (p^,fc(j,)) < ;^dist(y, J), 

where lq (^p^ ^y-^^ denotes the Sq-\ength of the curve p^ and C4 > is a constant, 

independent of Lq, depending only on the (asymptotic) angles between the curves 
p -r, ^ and J. 

By further shrinking the exceptional manifold J a little bit, and by selecting a 
suitably thin, one-sided neighborhood Ui = Us^ of J, we can achieve that the open 
25i-neighborhood of f/i (on the same side of J as Uq and Ui) is a subset of Uq- 
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For a varying 5, < S < Si, we introduce the layer 

Us = \ y e {Us\ Us/2) r\Do\3 a. sequence t„ 00 
(4.23) ^ 

such that S^'^y e {Us \ Us/2) for all n L 

Since almost every point of the layer {Us \ Us/2) HDq returns infinitely often to this 
set and the asymptotic equation 

f^{{Us\Us/2)nDo)^^-f^i{J) 
holds true, we get the asymptotic equation 



(4.24) fi{Us)^^-f^i{J). 



We will need the following subsets of Us- 



Us{oo) = {y e Us\ ki{y) = 00} . 

Here C3 is the constant from Corollary 4.20, the exact value of which will be specified 
later, at the end of the proof of Main Lemma 4.5. By selecting the pair of sets {Ui, J) 
small enough, we can assume that 



(4.26) ztubiy) > c^Si yyeUi. 

This inequality guarantees that the collision time %^(y)(?/) (y ^ Us) cannot be near 
any return time of y to the layer {Us \ Us/2), for 6 < di, provided that y e Usics). 
More precisely, the whole orbit segment S^~'^^~^^''^^^^^z will be disjoint from Ui, 
where z = S'y, t^^^y^_^{y) < t < H^^y^iy). 

Lemma 4.27. iJ,(Us\ Usic^)) = o{5) (small ordo of 5), as 5 ^ 0. 

Proof. The points y of the set Us \ Usics) have the property %^(y)_i(y) < C3. By 
doing another slight shrinking to J, the same way as in Note 4.18, we can achieve 
that %^(j,)(y) < 2c3 for all y E Us\ Usics), < d < Si. This means that all points 
of the set Us\Us{c3) are at most at the distance of S from the singularity set 
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U s-*{sn-). 

0<t<2c3 

This singularity set is a compact collection of codimension-one, smooth submani- 
folds (with boundaries), each of which is uniformly transversal to the manifold J. 
This uniform transversality follows from Lemma 4.3 above, and from the fact that 
the inverse images S~^{S7l~) {t > 0) of singularities can be smoothly foliated with 
local, concave orthogonal manifolds. Thus, the ^-neighborhood of this singularity 
set inside U s clearly has //-measure small ordo of 5, actually, of order < const -(J^. □ 

Lemma 4.28. For any point y e Us{oo) the curves Py,k{s) (0 < s < h{y, k)) have 
a C^-limiting curve Py^oois) (0 < s < h{y, oo)), with h{y, k) — > h{y, oo), as /c — > oo. 

Proof. Besides the concave, local orthogonal manifolds Sq(|/) = Sg'^d/) of (4.12) 
(where tl — t\{y) — \ {tk-i{y) + tfc(y))), let us also consider another type of con- 
cave, local orthogonal manifolds defined by the formula 

(4.29) So'(y) = S?(y) = SCy (5-*^ (^SCy^, {y' e M| q{y') = g(|/t-)})) , 

the so called "candle manifolds", containing the phase point y G Us{oo) in their 
interior. It was proved in §3 of [Cli(1982)] that the second fundamental forms 
B (Eo(j/), y) < are monotone non-increasing in k, while the second fundamental 

forms B ^Eq(|/), y^ < are monotone increasing in k, so that 

5(sg(y),y) <S(Eg(y),y) 
is always true. It is also proved in §3 of [Ch(1982)] that 

hm B (Sg(y), y) = hm B (Eg(y), y) := B^ < 

uniformly in y, and these two-sided, monotone curvature limits give rise to uniform 
C^-convergences 

lim Eg(y) = SS°(y), hm Sg(y) = SS°(y), 

t — >CXD t — >00 

and the limiting manifold T,Q^{y) is the local stable invariant manifold 'y^^\y) of y, 
once it contains y in its smooth part. These monotone, two-sided limit relations, 
together with the definition of the curves Py,ti — Py,k prove the existence of the 
C^-limiting curve Py,oo = hmfe_»oo Pj/,fc) h{y,k) — > h{y, oo), as A; — > oo. They also 
prove the inclusion Py,oo ([0, h{y, oo)]) C 7^*^(y). □ 
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Lemma 4.30. For any point y e Us{oo) the projection 11(2/) '■— Py,oo{h{ii,oo)) 
(e J) is either an element of the set A (of 4.4) or it is a forward singular point of J. 

Proof. Assume that the forward orbit of 11(2/) is non-singular. Select a return time 
t>cz for which S^y e Us C (Us \ U5/2) n Dq. The distance dist(5'*2/, J) between 
S^y and J is bigger than 5/2. According to the contraction result 4.20, if the 
contraction factor Lq^ is chosen small enough, the distance between (I\.{y)) and 
J stays bigger than 5/4, so (n(j/)) & Uq\J will be true. This means, however, 
that the forward orbit of 11(2/) is sufficient, according to (4) / (ii) of Definition 4.1. □ 

Corollary 4.31. ^{Us{oo)) =0. 

Proof. By our indirect assumption A = 0, so the previous lemma says that 

n(2/) G U {sn-) 

t>0 

for all y e Us{oo). Due to the transversality result 4.3 above, the singularity set 
Ut>o'^~* (»S7?.~) is a countable union of smooth, proper submanifolds of J. Since 
the curves Py, oo belong to (at least) the smoothness class, we get that Us{oo) 
is a countable union of proper, C^-smooth submanifolds of Us, hence fx (Us{oo)) = 
0. □ 

Next we need a useful upper estimation for the /U-measure of the set Usics) \ 
Us{oo) as 5 ^ 0. We will classify the points y G Us{cs) \ t/'j(oo) according to 
whether S^y returns to the layer [Us \ Us/2) H Dq (after first leaving it, of course) 
before the time t-j^_^^y^^_-^{y) or not. Thus, we define the sets 

Esic^) ={2/ G 1/5(03) \t^5(oo)| 3 < ti < t2 < %^(y)_i(?/) 

(4.32) such that S'^y ^Uo, S'^ye {U5 \ Us/2) n L>o}, 

Fsics) =Usic3) \ [Usioo) U Es{cs)] . 

Recall that the threshold t-j:^^y^_-^^{y), being a collision time, is far from any possible 

return time t2 to the layer (Us \ Us/2) H Do-, see the remark right after (4.26). 

Now we will be doing the "slight shrinking" trick of Note 4.18 the third (and 
last) time. We slightly further decrease J to obtain a smaller Ji with almost the 
same /Ui-measure. Indeed, by using property (6) of 4.1, inside the set J H B we 
choose a compact set Ki for which 

and no point of Ki ever returns to J. For each point x G Ki the distance between 
the orbit segment S^"'°''^^^x and J is at least e{x) > 0. Here is needed to guarantee 
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that we certainly drop the initial part of the orbit, which still stays near J, and 
C3 was chosen earlier. By the non-singularity of the orbit segment S^^°''^^^x and by 
continuity, the point x G Ki has an open ball neighborhood B{x) of radius r{x) > 
such that for every y G B{x) the orbit segment S^'^"''^^^y is non-singular and stays 
away from J by at least e{x)/2. Choose a finite covering UILi ^i^i) ^ -^1 -^i' 
replace J and Us hj their intersections with the above union (the same way as it 
was done in Note 4.18), and fix threshold value of 5o so that 

5o < ^ min{e(a;i) | z = 1, 2, . . . , n}. 

In the future we again keep the old notations J and Ug for these intersections. 
In this way we achieve that the following statement is true: 

Any return time t2 of any point y G (Us \ Us/2) H Do to 
(Us \ Us/2) n -Do is always greater than C3 for < 5 < 5i. 

For any phase point y G Esic^) we define the first return time t2 = ^2(2/) as the 
infimum of all the return times t2 of y featuring (4.32). By using this definition of 
t2{y), formulas (4.32)-(4.33), and the contraction result 4.20, we easily get 

Lemma 4.34. If the contraction coefficient Lq^ in 4.20 is chosen suitably small, 
then for any point y G £"5(03) the projected point 




(4.35) U{y):=p^-,^^^^{h{yMy)))eJ 

is either a forward singular point of J, or Il{y) belongs to the set A of regular 
sufficient points of J, defined in 4.4. 

Proof. Since t2{y) < t-^^(^y^_i{y), we get that, indeed, n(y) G J. Assume that the 
forward orbit of 11 (y) is non-singular. 

Since e {Us \ Us/2) n Dq, we get that dist (s^^^y^y, j) > 5/2. On the 

other hand, by using (4.33) and Corollary 4.20, we get that for a small enough 
contraction coefficient L^^ the distance between S^'^^^^y and ^'^^(f) (11 (j/)) is less 
than 5/4. (The argument is the same as in the proof of Lemma 4.30.) In this way 
we obtain that ^^^(s') (n(y)) eUo\J, so Il{y) G A, according to condition (4)/(u) 
in 4.1. □ 



Corollary 4.36. i^{Es{c3)) = 0. 
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Proof. The proof will be analogous with the proof of Corollary 4.31 above. Indeed, 
A = (according to our indirect assumption), and the exceptional set of Lemma 
4.2 (i. e. the set of all forward singular points of J) is a countable union of smooth, 
proper submanifolds of J. Since the curves Pyj^{y) belong to (at least) the 
smoothness class, we get that Es{cz) is a countable union of proper, C^-smooth 
submanifolds of Us, hence ^ {Es{cs)) = 0. □ 

For the points y e Fs{c3) we define the projection 11 (y) by the formula 

(4.37) U{y) := e dM. 

Now we prove 

Lemma 4.38. For the measure v (H (^5(03))) of the projected set H (^5(03)) C dM. 
we have the upper estimate 

Z/(n(F5(c3))) < 

where C2 > is the geometric constant (also denoted by C2) in Lemma 2 of [S- 
Ch(1987)] or in Lemma 4.10 of [K-S-Sz(1990)-I], C4 is the constant in (4.22) above, 
and 1/ is the natural T-invariant measure on that can be obtained by projecting 
the Liouville measure fi onto dM. along the billiard flow. 

Proof. Let y e Fs{c3). From the inequality t^_^(^y^_i{y) > C3 and from Corollary 

4.21 we conclude that Ztub (n(?/)) < C4Lq^S. This inequality, along with the fun- 
damental measure estimate of Lemma 2 of [S-Ch(1987)] (see also Lemma 4.10 in 
[K-S-Sz(1990)-I]) yield the required upper estimate for u {U (^5(03))). □ 

The next lemma claims that the projection H : Fs{cs) ^ dM (considered here 
only on the set ^5(03)) is "essentially one-to-one", from the point of view of the 
Poincare section. 

Lemma 4.39. Suppose that yi, 1/2 £ -^^(ca) are non-periodic points (5 < di), and 
n(yi) = n(y2)- We claim that yi and y2 belong to an orbit segment 5" of the billiard 
flow lying entirely in the one-sided neighborhood Uq — Us^ of J and, consequently, 
the length of the segment S is at most l.ldiam(J). 

Remark. We note that, obviously, in the length estimate l.ldiam(J) above, the 
coefficient 1.1 could be replaced by any number bigger than 1, provided that the 
parameter 5 > is small enough. 

Proof. The relation n(yi) = Il{y2) implies that yi and 2/2 belong to the same orbit, 
so we can assume, for example, that j/2 = S°'yi with some a > 0. We need to prove 
that S^^'"'^yi C Uq. Assume the opposite, i. e. that there is a number ti, < ti < a, 
such that S^^yi ^ Uq. This, and the relation S°-yi G (Us \ Us/2) nDo mean that the 
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first return of yi to (Us \ Us 72) nL>o occurs not later than at time t = a. On the other 
hand, since n(yi) = n(S'"yi) and yi is non-periodic, we get that > a, 

see (4.37). The obtained inequahty t-^_^(^y_^^_^{yi) > a > t2{y)i however, contradicts 
to the definition of the set ^5(03), to which yi belongs as an element, see (4.32). 
The upper estimate l.ldiam(J) for the length of S is an immediate corollary of the 
containment S G Uq. □ 

As a direct consequence of lemmas 4.38 and 4.39, we obtain 

Corollary 4.40. For all small enough 5 > the inequality 

l^{Fs{cs)) < I.IC2C4LQ ^5diam(J) 

holds true. 

Finishing the Indirect Proof of Main Lemma 4.5. 

It follows immediately from Lemma 4.27 and corollaries 4.31, 4.36, and 4.40 that 

fJ' {Us) < 1.2c2C4diam( J)Lq 

for all small enough S > 0. This fact, however, contradicts to (4.24) if Lq ^ is 
selected so small that 

1.2c2C4diam(J)LQ ^ < i/ii(J*), 

where J* stands for the original exceptional manifold before the three slight shrink- 
ings in the style of Note 4.18. Clearly, /xi(J) > (1 - 10-^)/xi( J*). The obtained 
contradiction finishes the indirect proof of Main Lemma 4.5. □ 



§5. The Theorem on Local Ergodicity Revisited 
Full Hyperbolicity Offsets the Ansatz 

In this section we are going to revisit the proof of the Theorem on Local Er- 
godicity (first proved by Chernov and Sinai as Theorem 5 in [S-Ch(1987)], see also 
Theorem 3.6 in [K-S-Sz(1990)-I] and Theorem 4.4 in [B-Ch-Sz-T(2002)]) in order 
to find a suitable, relaxed version of the Ansatz (see §2) sufficient for the proof of 
the mentioned theorem. This relaxed form of the Ansatz will be implied by the full 
hyperbolicity of the considered semi-dispersive billiard system. 

Assume, therefore, that the semi-dispersive billiard flow (M, {S'^jteR, //) is fully 
hyperbolic, i. e. all relevant Lyapunov exponents are nonzero ^u-almost everywhere. 
We want to prove the Theorem on Local ergodicity for (M, {S'^jigK, /i), that is, 
Theorem 5 of [S-Ch(1987)] without assuming the Ansatz, rather by postulating the 
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full hyperbolicity only. We will be constantly making references to the proof of 
Theorem 3.6 of [K-S-Sz(1990)-I]. 

We notice, first of all, that the generalization of Pesin's theory [P(1977)] is feasible 
for semi-dispersive billiards, see [K-S(1986)]. This theory asserts, among other 
things, that in a fully hyperbolic, semi-dispersive billiard flow (M, {^^jtgM, Z^) 

(1) the ergodic components have positive measure (hence there are at most count- 
ably many of them), and 

(2) the restriction of the flow (M, {S'^jteK, A*) to any of its ergodic components 
is K-mixing. 

Let is see now the part of the proof of Theorem 3.6 (of [K-S-Sz(1990)-I]) that 
actually uses the Chernov-Sinai Ansatz. It is the proof of the so called Tail Bound 
estimation, i. e. Lemma 6.1 there. We will be using the notations of that proof 
along with notations and definitions from the previous sections of this paper without 
risking to cause any conflict between these, slightly different systems of symbols. 
The only exceptions to this are that in the present section denotes the projection 
of the Liouville measure onto along the flow, just as in [K-S-Sz(1990)-I], and 
the constant C3 here is the same as the constant C3 of the cited paper, not our 
constant C3 in the previous section. 

In the current set-up, the exceptional manifold J is a smooth, round shaped piece 
of a post-singularity manifold. For technical reasons, the forward orbits S^^'^^^x of 
phase points x ^ Us will correspond to the backward orbits S^~^'^^ {S~'^ (T'^y)) 
{y G ^, notations of [K-S-Sz(1990)-I]), where x = S'^{—T'^y), and r runs over a 
fixed, suitably chosen, open time interval (a, h) with < a < 6. 

Note. It is clear that, by choosing the shape of the considered, small piece of post- 
singularity manifold J suitably, all postulates of 4.1 can be achieved to be true with, 
possibly, the only exception of (4) / (ii) . This property, however, will not be used in 
this section. We also note that (6) of 4.1 automatically holds true in this case since, 
according to Proposition 4.2, the set W has full measure in J. Finally, (5) of 4.1 is 
a direct consequence of (7.11) in [Sim(2003)]. 

Let us focus on the sets = U^^{5) and the corresponding numbers a^nm — 

ji {U^^^ from the proof of Lemma 6.2 of [K-S-Sz(1990)-I]. Let the ergodic com- 
ponents of the fiow (M, {^^jteiK, /i) be {Cfc| < A; < /cq}) where 1 < A; < 00. We 
remind the reader that, according to Lemma 4.10 of [K-S-Sz(1990)-I] and the Borel- 
Cantelli lemma, for /ii -almost every phase point x G 9M there exists a number 
C6 = cq{x) > such that ztubiT'^x) > C6(x)/n^ for all n G N or, equivalently for the 
flow, for /x-almost every phase point x G M there exists a number = cq{x) > 
such that Ztub{S^x) > CQ{x)/t'^ for alH > 1. This equivalence easily follows from 
the standard procedure of introducing the so called "transparent walls" for the flow, 
to make the map T have a finite horizon. We will be dealing with the flow and, to 
be definite, by cq = cq{x) we will mean the supremum (maximum) of all the above 
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constants. The key observation is that the sets U!^^^ — U^ ,^{6) can be spht fur- 
ther into countably many pieces, without hurting the apphcabihty of the arithmetic 
lemma 6.2, in the following way: 



Un,m,k,l,p = U^,m,k,l,p{^) ^ {v & U\ ZtubiT^'y) < Kn,cs5{y) ^C^^ 

(5.1) A"^ < «n,c35(y) < A"^+\ -T^y e Cfc, ^ < r(-T-y) < ^, 
^<C6(-T"y)<-}, 

for n, m, /, p G N, < /c < /cq, where r{—T'^y) is the radius of the largest open, 
geodesic ball lying fully inside the local stable manifold j^{—T"'y) and centered at 
—T'^y. (The so called inner geodetic radius of the local stable manifold 7*(— T"?/).) 

Note. In the case / = 1 or p = 1 we do not impose the unnecessary upper bound 
1/Z = 1 on r(-T"y) or 1/p = 1 on C6(-T"y). 

After this splitting of the basic sets C/^, the array of numbers 

J _ r-l (Tjh \ 

i,m,k,l,p "1 \ n,m,k,l,p) 



On, 



will be a 5-dimensional array, and the properties (i) and (iii) of Lemma 6.2 will 
obviously hold true: 

(i) J2n^n,m,k,l,p ^ ^m,k,l,p-, 
(iii) Em,k,l,pA^,k,l,p < oo. 

Here the numbers A^^^ i^p are independent of 5. In order to enable Lemma 6.2 
of [K-S-Sz(1990)-I] for aplication, the only thing we need to prove is the hypothesis 
(ii) of that lemma, i. e. 

Proposition 5.2, Relaxed form of the Ansatz. For every quadruplet (m, /c, l,p) 
{m,l,p eN, < k < ko) 

lim On m fc / o = 0- 

n>N 

We want to translate Proposition 5.2 to our language oi J, Us-, x — S'^{—T^y) G 
Us {a < T <b), and the language of the billiard flow via the correspondence 



(5.3) 



Vn,m,k,l,p = V^,m,k,l,P = '^^"'^^ (~^'' (^n,m,fc,Z,p) ) ^ ^s, 
bn,m,k,l,p = ^n,m,k,l,p ~ ^ {^n,m,k,l,p) — '^n,m,k,l,pi 
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for < 5 < 5i. For each ergodic component Ck select and fix an open region (in the 
relative topology of C/e) 7^ -Rfc C Cfc of sufficiency (hyperbolicity) in such a way 
that for every pair (m, A;) (m e N, < A; < ko) there exists a threshold N = N{m, k) 
such that 



(5.4) nn,cssix) > 

whenever n > N{m, k), x E Us C) Ck, and there is a time t, < t < t^n/2]{x)j such 
that S\Y{x)) n i?fc 7^ 0, where 

i^n,csd{x) is defined by 



(5-5) /^n, 035(2;) — I^n,cz5{]j) 1 

so that X = S'^{—T^y) for some r e (a, 6). Corresponding to (5.4) and the above 
said, we define the sets 

(5.6) 

for all t with <t < t\r,/2]ix), < r(x) < -, and < cq(x) < -\. 

In/^JV ^ ' - V p+ 1 ^ ' ~ p) 

Clearly, the monotonicity relations W^^^,^,;,^ 3 ^n2,k,i,p < na, and Vt^^;^,,,^ C 

^n,k,l,p for < hold true. 
Now we can prove 



Lemma 5.7. For every quintuplet {n,m,k,l,p) we have that V^,^^f^^i p C W^^^/j^/^p, 
whenever n > N{m, k). 

Proof. For any x e V;^,^ ^ we have 

where x = S'^{—T'^y) for some y E U and r e (a, 6). By putting this together with 
(5.4) and the assumption n > N{m,k), we get that S^{Y{x)) D Rk = $ for all t 
with < t < t[n/2]{x), i. e. X E W^^ a;,/,^ by the definition (5.6). □ 

Since the post-singularity set 

a<t<b 

can be covered by finitely many patches of the above sets J, instead of Proposition 
5.2 it is enough to prove 
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Proposition 5.8. For every quadruplet {m,k,l,p) 

1™ / bfi m k I V = ^■ 
5^0,Ar^oo ^ n,m,K,i,p 

n>N 

Proof of 5.8. We define the numbers 

(5.9) Cn,k,l,p = ci^k,l,p = V {Wn,k,l,p) ■ 

By the disjointness relations Vnj^^m,k,i,p^Vn2,m,k,i,p — ^or ni ^ n2 (see Lemma 6.3 
in [K-S-Sz(1990)-I]) and by the previous lemma we obtain that 

(5.10) Yl ^n,m,k,i,p < c%,k,i,p for N > N{m, k). 

n>N 

Thus, in order to prove 5.8 it is enough to prove 
Proposition 5.11. For every triplet {k,l,p) 

lim c%ui^ = 0. 

Proof of 5.11. Assume the opposite, i. e. that there exist two sequences 5i \ 0, 
Ni y oo, and a number a > such that 

(5.12) 4 > a for aU i e N. 
Define the sets 

(5.13) Gi = [j{Y{x)\xeW'^^^,^,^^}cCk. 

(We will suppress the dependence of the sets Gi on the fixed indices k, I, and p.) 
Observe that the sets Gi form a nested sequence 

(5.14) GiDG2DGsD ... 
and, by (5.12), 

(5.15) 5- V {Gi n Us, ) > a for z < j. 

An important point here is that the lower density estimation (5.15) can be extended 
to the case j < i (up to a constant multiplier), as the following lemma claims: 
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Lemma 5.16. There is a constant K = Kp > 1 (depending only on p) such that 

{Gi n Usj) > a/K for all i, j e N. 

Proof of 5.16. The local stable invariant manifolds 7^(0;) {x E Uq) are known to 
be uniformly transversal to the manifolds 

(5.17) Js = {x e dUs\ d(x, J) = 5} 

for 5 < 5o, see Lemma 4.3 above. We select and fix a smooth foliation 

(5.18) Uo = Us, = \J {Fs,,\ 0<S<So, ze K'^-^} , 

complementer and uniformly transversal to the foliation by the local stable mani- 
folds 7*( . ) with the additional property that 

(5.19) [j{Fs,,\zeR''-^}^Js 
ior < 5 < Sq. We remind the reader that 

d = dimFs^z = dimQ = dim7^(a;) + 1 = dimM — dim7*(a;) 

for X e Uq, and the pairs (5, z) form a, {d — l)-dimensional parameter domain. 

Now we take advantage of the fact that each constituent local stable manifold 
7*(.) of Gi contains a point x e k l p with < C6{x) < ^, hence the 
inequalities 

hold true for all points y G 7^(0;) with a constant K' > 1. It is a common feature 
of all proofs of the absolute continuity of the local stable foliation 7*( . ) (see, for 
example, Part II in [K-S(1986)]) that uniform lower and upper bounds are provided 
for the Jacobian of the holonomy map, once estimations of (5.20) are established. 
Therefore, there is a constant K — Kp > 1, depending on p and the constant K' of 

(5.20) , such that 

(5.21) K-^ < \\JacHs,z,S',z'{x)\\ < Kp 

for every point x e Fs^z ^ Gi. Here Hs^z,S',z' '■ F^^z ~^ Fg/^z' is the holonomy map 
defined by {H5^z,S',z'{x)} = 7*(a^) fl Fgi^z'- We remind the reader that, due to the 

'1' 2 ' 

size condition — — — < r(x) < y for all x e Gi, the holonomy map Hs^z,5',z'{ ■ ) is 

well defined for all points x G Fg^z H Gi for small enough 6 and 6'. Now (5.15) and 
(5.21) finish the proof of Lemma 5.16. □ 
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Finishing the proof of Proposition 5.11. 

Take the intersection Goo = fli^i of the nested sequence (5.14). By Lemma 
5.16 we have that 

57V (Goo nc/^,) > ^ 
j\p 

for all J, in particular n{Goo) > 0. For every constituent manifold 7*(y) C Goo 
and every i E N there exists a point Xi G Yiv) such that 5'*a;i -R^ for all t with 
< t < t[i/2]{x). Therefore, any point y e Goo will eventually avoid a shrunk 
ball inside Rk- This, however, contradicts to the ergodicity of the restricted system 
(Gfe, {<S'^}, A*|C'fe), thus finishing the proof of 5.11. □ 

The successful completion of the proof of 5.11 implies that the proof of the 
Theorem on Local Ergodicity works in every fully hyperbolic semi-dispersive billiard 
system without postulating the Chernov- Sinai Ansatz. □ 

An immediate consequence of the results of this section is 

Corollary 5.22. If a semi-dispersive billiard system is fully hyperbolic, then all of 

its crgodic components are open (modulo the zero sets, of course). In particular, 
by virtue of the Theorem in [Sim(2002)] (see the introduction of that paper), every 
hard ball system has open ergodic components. 



6. Proof of Ergodicity 
The Grand Induction 



By using several results of Sinai [Sin(1970)], Chernov-Sinai [S-Ch(1987)], and 
Kramli-Simanyi-Szasz, in this section we finally prove the ergodicity (hence also 
the Bernoulli property, see Chernov-Haskell [C-H(1996)] or Ornstein- Weiss [O- 
W(1998)]) for every hard ball system (M, {/S*}, //) by carrying out an induction 
on the number N (> 2) of interacting balls. But before carrying out this induction, 
let us ask the naturally arising question: Why do we need the induction hypoth- 
esis at all? Why cannot we directly apply our Corollary 5.22 above to conclude 
that (M, {S^}, jj.) has open ergodic components, and than disprove the possibility 
of more than one component by directly using Main Lemma 4.5? The answer is 
this: The applicability of 4.5 requires that the set separating distinct ergodic com- 
ponents be a smooth exceptional manifold of non-sufficiency! It simply does not 
work if highly complicated, closed sets of eventually splitting orbits also participate 
in the business of separating different ergodic components. That is why Theorem 
5.1 of [Sim(1992)-I] needs to be used in order to guarantee that the set of even- 
tually splitting orbits is slim, hence incapable of separating distinct open ergodic 
components. The cited theorem, however, uses the induction hypothesis. 
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The base of the induction (i. e. the ergodicity of any two-ball system on a flat 
torus) was proved in [Sin(1970)] and [S-Ch(1987)]. 

Assume now that (M, {S^}, fi) is a given system of A'^ (^3) hard spheres with 
masses mi, m2, . . . , mjv and radius r > on the fiat unit torus = M'^/Z'^ (z/ > 2), 
as defined in §2. Assume further that the ergodicity of every such system is already 
proved to be true for any number of balls N' with 2 < N' < N. We will carry out 
the induction step by following the strategy for the proof laid down in the series of 
papers [K-S-Sz(1989)], [K-S-Sz(1990)-I], [K-S-Sz(1991)], and [K-S-Sz(1992)]. 

By using the induction hypothesis, Theorem 5.1 of [Sim(1992)-I], together with 
the slimness of the set A2 of doubly singular phase points, shows that there exists 
a slim subset S'l C M of the phase space such that for every x G M \ S'l the point 
X has at most one singularity on its entire orbit S^~°^'°°^x, and each branch of 
X is not eventually splitting in any of the time directions. By Corollary 
3.26 and Lemma 4.2 of [Sim(2002)] there exists a locally finite (hence countable) 
family of codimension-one, smooth, exceptional submanifolds Jj C M such that for 
every point x ^ (Ui <^») ^ "^i ^^e orbit of x is sufficient (geometrically hyperbolic). 
This means, in particular, that the considered hard ball system (M, {5'*}tg]R, fi) is 
fully hyperbolic. According to our result in §5, the Theorem on Local Ergodicity 
by Chernov and Sinai (Theorem 5 in [S-Ch(1987)], see also CoroUary 3.12 in [K-S- 
Sz(1990)-I]) is applicable at any hyperbolic (sufficient) phase point x E M. with at 
most one (simple) singularity on its orbit. By that theorem and the main result of 
[B-Ch-Sz-T(2002)]), an open neighborhood 1/^ 3 x of any phase point x ^ (|J^ Ji) U 
belongs to a single ergodic component of the billiard fiow. (Modulo the zero 
sets, of course.) Therefore, the billiard flow {S^} has at most countably many, open 
ergodic components Ci, C2, .... 

Assume that, contrary to the statement of our theorem, the number of ergodic 
components Ci, C2, ... is more than one. The above argument shows that, in 
this case, there exists a codimension-one, smooth (actually analytic) submanifold 
J C M \ dM. separating two different ergodic components Ci and C2, lying on the 
two sides of J. By the Theorem on Local Ergodicity for semi-dispersive billiards, no 
point of J has a sufficient orbit. (Recall that sufficiency is clearly an open property, 
so the existence of a sufficient point y E J would imply the existence of a sufficient 
point y' E J with a non-singular orbit.) By shrinking J, if necessary, we can achieve 
that the inflnitesimal Lyapunov function Q{n) be separated from zero on J, where 
n is a unit normal fleld of J. By replacing J with its time-reversed copy 

-J = {(g,^)eM| (g,-^)e J}, 

if necessary, we can always achieve that Q{n) < —c\ < uniformly across J. 

There could be, however, a little difficulty in achieving the inequality Q{n) < 
across J. Namely, it may happen that Q{nt) = for every t e R. This is, however, 
shown to be impossible in Remark 7.9 of [Sim(2003)]. 
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To make sure that the submanifold J is neatly shaped (i. e. it fulfills (2) of 4.1) 
is a triviality. Condition (3) of 4.1 clearly holds true. We can achieve (4) as follows: 
Select a base point xq & J with a non-singular and not eventually splitting forward 
orbit S^^'°^^xo. This can be done according to the transversality result Lemma 
4.3 (see also 7.12 in [Sim(2003)]), and by using the fact that the points with an 
eventually splitting forward orbit form a slim set in M (Theorem 5.1 of [Sim(1992)- 
I]), henceforth a set of first category in J. After this, choose a large enough time 
T > so that S'^xq ^ 9M, and the symbolic collision sequence Eq = E (S'^'^'^^xq) 
is combinatorially rich in the sense of Definition 3.28 of [Sim(2002)]. By further 
shrinking J, if necessary, we can assume that S^{J) n dM = and is smooth 
on J. Choose a thin, tubular neighborhood Uq of J in M in such a way that S'^ be 
still smooth across Uq, and define the set 



(6.1) NS (Uo, So) = e C/o| dimWo (^S^^'^^x^ > l} 

of not Eo-sufficient phase points in Uq. Clearly, J C NS (^o, Sq). We can assume 

that the selected (generic) base point Xq E J belongs to the smooth part of the 

closed algebraic set NS (Uq, Sq). This guarantees that actually J = NS (Uq^ Eq), 

as long as the manifold J and its tubular neighborhood are selected small enough, 
thus achieving property (4) of 4.1. 

Proof for why property (6) of Definition 4.1 can be assumed. 

We recall that J is a codimension-one, smooth manifold of non-sufficient phase 
points separating two open ergodic components, as described in (0)-(3) at the end 
of §3 of [Sim(2003)]. 

Let P be the subset of J containing all points with non-singular forward orbit 
and recurring to J infinitely many times. 

Lemma 6.2. iJ,i{P) = 0. 

Proof. Assume that Hi{P) > 0. Take a suitable Poincare section to make the 
time discrete, and consider the first return map T : P — > P of P. Let Pn be the 
set of all points x G P for which the time it takes to return to P belongs to the 
interval [n, n + 1), n = 0, 1, 2,. ... The sets Pn provide a measurable partition 
of P. Observe that the measurable subsets T{Pn) of P are mutually disjoint, and 
/ii(T(Pn)) > idi{Pn) (for all n for which HiiPn) > 0) by the hypersurface measure 
expansion theorem proved in [Ch-S(2006)]. The obtained contradiction proves the 
lemma. □ 

Next, we claim that the above lemma is enough for our purposes to prove (6) of 
4.1. Indeed, the set W G J consisting of all points x E J never again returning to J 
(after leaving it first, of course) has positive //i-measure by the Lemma 6.2. Select 
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a Lebesgue density base point xq & W for W with a non-singular forward orbit, 
and shrink J at the very beginning to such a small size around xq that the relative 
measure of in J be bigger than 1 — 10~®. 

Finally, Main Lemma 4.5 asserts that A 7^ 0, contradicting to our earlier state- 
ment that no point of J is sufficient. The obtained contradiction completes the 
inductive step of the proof of the Theorem. □ 

Appendix. The Constants of §3-4 

In order to make the reading of sections 3-4 easier, here we briefly overview the 
hierarchy of the constants used in those sections. 

1. The geometric constant cq > of Proposition 3.5 is a lower estimation for the 
"curvature" {Sqo, 6vo)/\\6qo\\^ of an expanding tangent vector {Sqo, Svq). 

2. The geometric constant — ci < provides an upper estimation for the infin- 
itesimal Lyapunov function Q{n) of J in (5) of Definition 4.1. It cannot be freely 
chosen in the proof of Main Lemma 4.5. 

3. The constant C2 > is present in the upper measure estimation of Lemma 
2 of [S-Ch(1987)], or Lemma 4.10 in [K-S-Sz(1990)-I]. It cannot be changed in the 
course of the proof of Main Lemma 4.5. 

4. The contraction coefficient < Lq^ « 1 plays a role all over §4. It must 
be chosen suitably small by selecting the time threshold C3 >> 1 large enough (see 
Corollary 4.20), after having fixed Uq, Sq, and J. The phrase "suitably small" for 
Lq^ means that the inequality 

^ 0.25//i(J*) 
° 1.2c2C4diam(J) 

should be true, see the end of §4. 

5. The geometric constant C4 > of (4.22) bridges the gap between two distances: 
the distance dist(y, J) between a point y G Us and J, and the arc length Iq (^p^ ^y-^^- 
It cannot be freely chosen during the proof of Main Lemma 4.5. 
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